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1.0  INTRODUCTION 


The  usefulness  of  holography  as  a  data-gathering  tool  is  determined 
primarily  by  the  accuracy  with  which  data  can  be  reconstructed  from  the 
hologram,  A  study  of  this  problem  can  be  approached  by  applying  a  linear 
systems  analysis  to  the  holographic  process  (Refs.  1  and  2), 

By  the  principle  of  superposition  for  linear  systems,  the  presence 
of  any  one  excitation  is  not  affected  by  the  response  attributable  to  other 
excitations.  A  complicated  input  function  can  be  considered  to  be  com¬ 
posed  of  multiple  simpler  functions,  such  as  the  components  of  a  Taylor 
series  or  a  Fourier  spectrum,  and  the  response  of  the  system  can  be 
studied  for  each  component.  One  such  approach  is  to  study  the  transform 
characteristics  of  the  system  to  an  impulse  function  input.  The  output 
caused  by  a  complicated  input  is  then  the  super  positioning  integral 
(mathematically,  the  convolution)  of  the  impulse  transform  function  over 
the  input  function  domain. 

The  optical  analog  of  an  impxflse  function  is  a  point  source  of 
radiation.  Invoking  Huygen's  principle,  a  radiating  surface  illuminance 
is  the  integration  of  a  point  source  over  the  radiating  area.  By  consid¬ 
ering  the  transfer  response  of  an  optical  system  to  a  point  radiating 
source,  one  can  find  the  essential  features  of  the  resolution  character¬ 
istics  of  the  system. 

Point  sources  of  radiation  are  described  in  scalar  diffraction  theory 
by  Green's  function 


G 


(1) 


where  A  is  the  amplitude,  k  =  2  tt/  X  ,  X  is  the  wavelength,  and  r  is  the 
distance  between  the  point  source  and  the  point  of  observation  (i.  e.  , 
detection).  The  physical  analog  of  Green's  function  is  the  E  field  obtained 
in  the  far  field  radiation  zone  of  an  electric  dipole  (Ref.  3): 

E  =  k2  (S  ^  P  X  (2) 

where  p  is  the  dipole  moment  and  A  is  the  unit  vector.  This  further 
suggests  that  scalar  diffraction  theory  should  be  restricted  to  the  far 
field  radiation  zone  requirements.  In  other  words,  the  radiation  wave¬ 
length  should  be  small  compared  to  the  diffraction  aperture  which,  in 
turn,  should  be  small  compared  to  the  observation  distance. 
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It  is  necessary  to  point  out  that  the  Green's  function  description 
is  not  adequate  to  describe  plane  wave  radiation.  Plane  waves  would,  by 
Green's  function,  need  to  originate  from  a  point  infinitely  distant  from  the 
observation  point;  this  would  reduce  its  amplitude  to  zero.  A  plane  wave 
being  characterized  by  its  amplitude  and  phase  front  is  adequately 
described  by 


G,  -  B.W  (3) 

where  the  phase  W  is  defined  relative  to  some  plane  perpendicular  to  the 
propagating  direction  along  the  radiation  path. 

By  the  Poynting  theorem,  the  observed  intensity  is  obtained  from 
the  product  of  the  radiating  field  times  its  complex  conjugate.  By 
requiring  that  the  intensities  obtained  from  the  plane  wave  and  point  source 
radiation  have  the  same  units,  the  units  for  A  and  B  are,  respectively, 
determined. 

2.0  HOLOGRAPHIC  RECORDING  AND  RECONSTRUCTION 


The  parameters  for  recording  and  reconstructing  the  hologram  of 
a  point  source  are  defined  in  Fig.  1. 

yi 


Figure  1.  Definition  of  parameters. 
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The  (X|,  with  i  =  0,  1,2  are,  respectively,  the  point  source  plane,  the 
holograna  recording  plane,  and  the  reconstructed  image  plane.  Let 
and  X2  be  the  wavelength  of  the  light  in  the  recording  and  reconstruction 
mode,  and  let  s^  and  Sg  be  the  distance  between  the  point  source  plane 
and  the  hologram  plane  and  the  reconstructed  image  plane,  respectively. 

The  recording  is  made  by  combining  in  the  (xi,y]^)  plane  light 
from  the  point  source  located  at  (xq  =  0,  Jq  =  0)  and  plane  waves 
propagating  in  the  positive  z -direction.  After  suitably  developing  the 
holographic  film,  the  reconstruction  is  made  by  illuminating  the  hologram 
in  the  (x2,yj^)  plane  by  plane  waves  again  propagating  along  the  positive 
z-axis.  All  the  optical  axes  have  been  made  coincident  for  mathematical 
simplicity,  not  because  of  any  inherent  limitation  in  the  holographic 
process. 

Chosing  the  (xq,  Yq)  plane  for  zero  reference  phase  of  the  plane 
waves, then  W  =  0  and  the  field  U  in  the  hologram  plane  will  be 


U  =  B  + 


(4) 


Restricting  the  field  to  a  paraxial  region,  one  can  use  the  binominal 
expansion. 


(5) 


where  rf  =  +  Xi  .  The  assumption  is  made  that  higher  order  terms 

are  small  compared  to  so  that  they  contribute  negligibly  to  the  phase 
in  U.  The  intensity  in  the  hologram  plane  is  obtained  from  the  product 
of  the  field  times  its  complex  conjugate  which  gives 


B" 


1  + 


2  „2 


=  Co  +  Cj 


(6) 


This  result  is  similar  to  a  Fresnel  zone  plate  (Ref.  4).  The  circular 
symmetry  about  the  z-axis  is  important  in  the  reconstruction  process. 
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The  cosine  term  is  the  information  carrying  term;  the  Cq  term  contributes 
a  constant  bias  term  to  the  film  exposure  level  which  can  generally  be 
ignored  in  the  reconstruction. 

The  reconstructed  amplitude  field  Uj^(x2,y2)  resiiLting  from  an 
illuminated  hologram  field  Upj(x]^jy]^)  is 

Ur  =  - 

The  field  Ujj  is  the  product  of  an  illuminating  radiation  source  and  the 
normalized  hologram  transmittance.  When  random  phase  shifts  attrib¬ 
utable  to  random  film  emtilsion  thicknesses  can  be  neglected,  the  complex 
amplitude  transmittance  of  the  hologram  is  a  linear  mapping  of  the 
complex  amplitude  incident  during  exposure,  provided  the  film  exposure 
level  is  chosen  properly  (Ref.  2).  Then  for  plane  incident  waves  B' 
on  the  hologram. 


where  and  C  ^  have  absorbed  the  normalization  and  the  film  transmit¬ 
tance  characteristics.  With  the  same  paraxial  approximation  used  as 
before 


cos  (  n,  S2)  dxjdyj 


(7) 


sg  =  [zf  +  (xj  -  xg)^  +  (yj  -  y2)^]^* 

^  '’2  “  ^*‘l'‘2  ~ 

These  approximations  also  allow  using  cos(n,  s)  —  1.  Changing  to  polar 
coordinates  and  using  the  integral  definition  of  the  Bessel  function,  one 
obtains 


277  J>)  = 


(10) 


then  approximating  that  S2  "  for  the  denominator  of  Green's  function. 


U„  = 


'2^1 


ik2r2/2zi  /*‘2‘'2T\  **^2 (t^ ^ ’'l)ridr j 


(11) 
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R  is  the  maximum  aperture  determined  by  the  size  of  the  hologram  or 
by  the  amplitude  sensitivity  of  the  film  emulsion.  The  integration  limits 
can  be  made  constant  by  using  the  change  of  variable  r^^  =  Rw;  then 


/2zi 


2^1 


i(k2R/2zpw^ 


»dw 


(12) 


After  expanding  the  cosine  term  of  Ujj  the  reconstruction  amplitude  can 
be  written  as 


3  1  iu'W^ 

Uj{  =  2  K./  Jo{vw)  e  J  wdw 

j=l  ^0 

where 


and 


(13) 


(14) 


The  three  integrals  are  successively  the  reconstruction  caused  by  a 
d-c  bias  on  the  film,  the  virtual  image,  and  the  real  image.  Following 
from  the  recording  geometry  described  in  Fig.  1.  ,  all  three  images 
overlap.  The  spatial  separation  of  the  images  obtained  with  an  off- 
axis  (side  band)  holographic  system  would  change  the  details  but  not 
the  general  form  of  the  reconstruction  integral.  Therefore,  the  integral 
contains  the  salient  features  of  a  point  source  reconstruction. 
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3.0  RECONSTRUCTION  CHARACTERISTICS 


The  solution  of  the  integral  of  the  form 

V  =  J  Jj^(vw)  wdw  (15) 

0 

can  be  made  using  Lommel  functions  (Ref.  5).  The  intensity  obtained 
from  YY'^  can  be  graphically  represented  by  a  set  of  isophote  diagrams 
(Refs.  5  and  6)  in  which  contour  lines  of  equal  intensity  near  the  focus 
of  the  optical  system  are  plotted  on  a  (u,  v)  coordinate  system.  The 
importance  of  these  isophote  diagrams  is  that  the  Eq.  (15)  solution 
they  represent  was  worked  out  in  analyzing  the  three-dimensional  light 
distribution  near  the  focus  of  a  diffraction-limited  lens.  In  other  words, 
the  intensity  field  in  the  focused  region  of  the  reconstruction  of  the 
hologram  of  a  point  source  is  identically  the  same  as  the  intensity  field 
near  the  focal  volume  of  a  diffraction  limited  lens.  Thus,  the  focusing 
characteristics  of  the  lens  apply  equally  well  to  the  holographic 
reconstruction.  In  particular, 

1.  The  fraction  of  the.  total  energy  which  falls  within  a 
given  area  about  the  optical  axis  for  any  constant  u 

is  the  same  for  point  images  obtained  by  lens  projection 
and  holographic  reconstruction  (where  the  reconstructed 
images  are  separated); 

2.  The  phase  distribution  near  focus  is  the  same  in  both 
cases,  i.  e.  ,  in  the  focal  waist  the  phase  surfaces 
propagate  like  plane  waves. 

Two  special  solutions  to  Eq.  (15)  are  available.  When  the  focus 
condition  is  satisfied,  ^2^1  =  -^i^o  (the  +(-)  is  for  the  real  (virtual)  image) 
then  u  =  o  and 

1 

U  =  /  J^(vw)wdw  =  - -  (16) 

0  ^ 

The  resulting  intensity  (normalized  to  unity  at  the  focal  center)  becomes 


(17) 
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which  is  the  same  as  the  Fraunhofer  diffraction  from  a  circular  aperture 
of  radius  R.  Analogously,  the  hologram  behaves  like  a  lens  of  radius  R 
with  a  focal  distance  z. 


The  second  special  solution  is  obtained  when  r2  -  o  which  is 
satisfied  along  the  optical  axis.  Then  v  =  o  and 

and  the  normalized  intensity  becomes 


2 


Sine 


(19) 


where  sine  (a)  =  sin  (a)/ a  and  the  +(-)  represents  the  virtual  (real)  image. 
This  intensity  distribution  can  be  considered  as  a  focal  depth  spread 
function  of  a  point  image.  In  other  words,  it  gives  the  intensities  at 
points  a  small  distance,  Az,  away  from  the  true  focus  position.  Consider 


then  approximately 


+  A2Z1  + 


Sine 


n-R^  Az 


4A 


2^1 


(20) 


(21) 


The  out-of -focus  distance  is  Az/2.  The  focal  depth  can  be  considered  to 
be  the  distance  Az  from  out-of-focus  position  (-Az/2)  to  out-of-focus 
position  (+Az/2)i  such  that  the  intensity  does  not  drop  below  a  given  part 
p  of  the  maximum  (at  Az  =0).  Usually  a  20 -percent  reduction  from  the 
maximum  intensity  (i.  e.  ,  p  =  0.  2)  is  allowed  (Ref.  5).  By  expanding  the 
sine  function  and  using  the  first  two  terms,  it  follows  that 

=  ^2  (y)  (22) 


The  spread  between  the  80 -percent  intensity  points  on  either  side  of 
focus  is 


(23) 
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The  lens  analogy  is  that  the  focal  depth  of  the  lens  increases  with  the 
f-number,  represented  by  z^/R. 


4.0  FOCAL  DEPTH  EXPERIMENT 


The  effectiveness  of  Eq.  (23)  was  tested  experimentally.  An 
object  field  of  equally  spaced  100-p  -diam  wires  was  made.  The  object 
is  pictorially  like  a  ladder  in  which  the  wires  serve  as  rungs.  The 
wire  spacing  was  3.  048  ±  0.  008  cm.  The  wires  were  positioned  in  a 
5 -cm  collimated  He-Ne  laser  beam  such  that  the  longitudinal  axis  of 
the  "ladder"  was  displaced  from  the  illumination  propagation  direction 
sufficiently  that  all  the  wires  were  illuminated.  The  effective  spacing 
between  the  wires  along  the  direction  of  the  beam  propagation  was 
3.  043  ±  0.  008  cm.  Holograms  were  made  using  Agfa  10E75  film.  The 
distance  from  the  nearest  wire  relative  to  the  film  plane  was  32  cm. 
Reconstruction  was  by  a  collimated  He-Ne  laser  beam.  In  the  recon¬ 
struction  the  distances  between  successive  wire  images  (judged  to  be 
at  their  best  focus)  were  measured.  By  averaging  and  calculating  the 
standard  deviation,  the  distance  between  the  reconstructed  images  was 
found  to  be  3.  04  ±  0.  15  cm.  From  multiple  repetition  of  the  experiment, 
the  standard  deviation  was  found  to  vary  by  about  ±0,  03  cm.  In  corre¬ 
lating  the  separation  distance  with  the  uncertainty  of  the  focal  position 
of  each  reconstructed  wire  image  it  was  assumed  that  each  wire  con¬ 
tributes  an  equal  amount  to  the  uncertainty  of  the  distance  between 
the  wire  image.  Then  the  image  position  uncertainty  is  one  half  the 
uncertainty  of  the  separation  distance.  Thus,  the  focused  image  is 
located  with  an  uncertainty  of  one  half  the  above  standard  deviation; 
in  other  words,  Az  =  (0.  15  ±  0.  03)  cm.  It  should  be  emphasized 
that  this  is  an  average  from  many  data  points;  it  is  not  intended  to 
imply  that  each  wire  reconstruction  can  be  focused  satisfactorily 
within  such  a  tolerance.  Rather,  the  average  result  implies  that  any 
particular  focused  image  location  cannot  be  measured  with  greater 
accuracy. 

Equation  (23)  cannot  be  applied  directly  to  finding  the  theoretical 
image  position  spread  for  this  experiment  because  z  varies  for  each 
wire.  However,  analogously  to  the  experimental  averaging  procedure, 
one  can  find  an  average  Az  spread  due  to  all  the  wire  images.  When  R 
is  considered  a  constant,  the  averaging  procedure  is  facilitated  by 
letting 

Zj  =  z  +  md  (24) 
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The  term  z  is  the  distance  from  the  film  plane  to  the  nearest  wire  (32  cm), 
d  is  the  spacing  between  wires  (3.  043  cm),  and  m  ranges  from  zero  to  N 
the  total  number  of  wires  in  the  object  field.  An  average  Az  for  the  80% 
intensity  envelope  case  is  obtained  from 


N-l 


s 


m=o 


Az 


n 


I  ^ 

N  r2 


2  (z  +  nd)^ 


+  zd(N-l)  +  -d-^tN 
6 


1)(2N 


-»] 


(25) 


There  remains  the  problem  of  determining  an  effective  aperture 
radius  R  of  the  hologram.  This  parameter  is  usually  described  in  terms 
of  the  resolution  characteristics  of  the  film  used,  the  assumption  being 
that  oscillations  in  the  diffraction  pattern  are  recorded  up  to  the  limit  of 
the  number  of  lines  per  millimeter  which  the  film  can  resolve.  It  is  easy 
to  show  that  this  assumption  does  not  satisfactorily  relate  to  the  problem 
at  hand.  Consider  the  diffraction  from  a  single  wire.  The  oscillatory 
part  of  the  intensity  varies  as 


(26) 


where  a  is  the  wire  cross-sectional  radius  and  the  other  parameters  are 
as  defined  before.  The  sine  term  with  the  linear  x-y  argument  serves  as 
an  envelope  to  the  other  more  rapidly  oscillating  sine  term.  It  can  be 
shown  that  if  a  film  has  a  resolution  given  as  n  lines  per  millimeter, 
the  distance  in  centimeters  for  which  this  limit  is  reached  by  the  wire 
diffraction  is 


Xj  =  lOnAjz^  (27) 

For  typical  parameters:  =  0.  6328iu  (He-Ne  laser),  n  =  2000  (Agfa 

Gaevert  10E75),  Zq  =  5  cm;  then  Xj^  =  63.3  cm.  Obviously,  this  far 
exceeds  the  size  of  film  plates.  Furthermore,  the  paraxial  approxi¬ 
mations  are  so  thoroughly  violated  that  the  accuracies  of  Eqs.  (23  and 
26)  are  very  questionable. 

The  effective  aperture  radius,  R,  could  in  principle  be  determined 
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by  knowing  the  dynamic  range  of  the  film  and  then  calculating  the  smallest 
visibility  which  could  be  recorded.  This  information  was  not  available 
and  consequently  a  more  direct,  although  less  rigorous,  approach  was 
used.  The  intensities  of  the  diffraction  pattern  of  the  hologram  of  a 
single  wire  were  obtained  by  scanning  the  hologram  with  a  microdensitom¬ 
eter.  It  was  found  that  the  oscillatory  information  trace  could  not  be 
distinguished  from  random  variations  past  approximately  the  fourth 
minimum  of  the  linear  argument  sine  term.  This  corresponds  to  a 
distance  from  the  center  of  the  diffraction  pattern  of 


4^0 


(28) 


This  distance  is  equated  to  R.  It  may  appear  that  R  shoiald  be  equal  to 
this  is  not  the  case  because  it  may  be  shown  from  the  mathematics 
of  the  reconstruction  of  the  wire  hologram  that  the  diffraction  on  each 
side  of  the  center  of  the  pattern  contributes  only  to  the  reconstruction 
of  the  edge  of  the  respective  side  of  the  wire.  It  is  as  though  the  hologram 
were  made  of  two  holograms,  each  reconstructing  half  of  the  wire  image. 


When  the  reconstruction  is  in  focus  and 


R  = 


44 


(29) 


Substituting  Eq.  (29)  into  Eq,  (23),  z^  now  appears  in  both  the  numerator 
and  denominator,  therefore  cancelling,  and  the  averaging  processing  of 
Eq.  (25)  is  no  longer  necessary.  The  depth  of  field  becomes  simply 


(30) 


The  depth  of  focus  in  dependent  only  on  the  wire  radius  and  on  a  change 
of  wavelength  between  recording  and  reconstruction.  For  the  experi¬ 
mental  conditions  of  Xj  =  =  0.  6328  ju  and  a  =  50  M  ,  one  obtains 

Az  =  1.  0  mm.  This  is  satisfyingly  close  to  the  (1.  5  ±  0.  3)  mm  depth 
of  field  observed  experimentally. 


The  above  manner  of  obtaining  the  effective  film  aperture,  R, 
rests  on  judgmental  criterion  (i.  e.  ,  visual  interpretation  of  the  micro¬ 
densitometer  trace)  rather  than  more  definitive  characteristics  based 
on  film  response  and  development  effects.  In  that  sense  the  above 
result  should  be  only  cautiously  interpreted  as  a  successful  experimental 
verification  of  the  theory.  Further  measurements  of  film  dynamic 
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range  and  visibility  function  are  necessary  to  justify  more  firmly  the 
above  determination  of  effective  film  apertures. 

An  additional  interpretational  difficulty  of  the  experimental  results 
is  that  the  intensity  range  of  the  focal  depth  was  rather  arbitrarily  chosen 
for  20 -percent  variations  from  the  maximum  (in  focus)  value.  This  value 
is  given  by  Born  and  Wolf  (Ref.  5)  but  without  any  discussion.  Different 
values  will,  of  course,  change  the  theoretical  depth  of  focus  as  can  be 
seen  from  Eq.  (22). 


5.0  SUMMARY 


The  focused  image  of  the  reconstruction  of  a  hologram  of  a  point 
radiation  source  is  shown  to  have  identically  the  same  intensity  and  phase 
characteristics  as  a  point  image  obtained  by  a  diffraction  limited  lens. 
Specifically  the  reconstructed  image  has  an  intensity  spread  function 
(Eq.  17)  and  a  focal  depth  (from  Eq.  21),  which  are  the  same  as  for  a 
diffraction  limited  lens,  and  within  its  focal  waist  the  propagation  is 
plane.  That  the  resiilts  of  the  two  imaging  processes  yield  identical 
results  is  interesting  because  their  focusing  phenomena  are  different, 
one  being  a  refractive  process,  the  other  a  diffractive  process.  The 
similarity  of  the  results  becomes  less  striking  when  the  diffraction 
limited  lens  derivation  is  examined  in  its  philosophical  components. 

In  that  derivation  it  is  assumed  that  perfectly  converging  coherent 
spherical  waves  (produced  by  the  refractive  phenomenon)  are  limited  in 
extent  by  an  aperture.  The  subsequent  results  are  all  due  to  diffraction 
by  the  aperture.  Essentially  the  refractive  phenomenon  is  not  really 
an  integral  part  of  the  derivation.  Nevertheless,  its  results  have  been 
well  born  out  by  experimental  studies  (Ref.  7)  and  its  acceptance  is 
universal. 

An  experiment  was  performed  to  test  the  focal  depth  of  the 
reconstructed  image.  The  theoretical  predictions  were  observed 
experimentally  within  the  cautions  of  using  interpretational  data. 

These  resiilts  show  that  depth  of  focus  is  dependent  on  the  object  size. 

The  size  of  the  objects  chosen  (wires)  were  designed  to  simulate  typical 
particle  field.  To  this  effect,  it  is  shown  that  the  depth  of  focus  of 
particles  of  50 radius  cannot  be  better  than  about  one  millimeter. 
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